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Abstract 

Stochastic Loewner evolution (SLE) is a differential equation driven by a one-dimensional Brow- 
nian motion (BM), whose solution gives a stochastic process of conformal transformation on the 
upper half complex-plane EL As an evolutionary boundary of image of the transformation, a ran- 
dom curve (the SLE curve) is generated, which is starting from the origin and running in EI toward 
the infinity as time is going. The SLE curves provides a variety of statistical ensembles of impor- 
tant fractal curves, if we change the diffusion constant of the driving BM. In the present paper, 
we consider the Schwarz-Christoffel transformation (SCT), which is a conformal map from EI to 
the region EI with a slit starting from the origin. We prepare a binomial system of SCTs, one of 
which generates a slit in EI with an angle an from the positive direction of the real axis, and the 
other of which with an angle (1 — o)tt. One parameter k > is introduced to control the value of 
a and the length of slit. Driven by a one-dimensional random walk, which is a binomial stochastic 
process, a random iteration of SCTs is performed. By interpolating tips of slits by straight lines, 
we have a random path in H, which we call an Iterative SCT (ISCT) path. It is well-known that, 
as the number of steps N of random walk goes infinity, each path of random walk divided by 
y/N converges to a Brownian curve. Then we expect that the ISCT paths divided by (the 
rescaled ISCT paths) converge to the SLE curves in N — > oo. Our numerical study implies that, 
for sufficiently large N, the rescaled ISCT paths will have the same statistical properties as the 
SLE curves have, supporting our expectation. 

PACS numbers: 05.40.-a, 05.45.Df, 02.30.-f 
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I. INTRODUCTION 



One of the highest topics of recent progress in statistical physics of critical phenomena 
and random fractal patterns is introduction of the Stochastic Loewner Evolution (SLE) 

nn 

by Schramm [1H3J. The SLE will provide a unified theory of statistics of random curves 
in the plane, which covers (continuum limits of) random interfaces characterizing surface 
critical phenomena in equilibrium (e.g. the percolation exploration process, the critical Ising 
interface), models of random chains in polymer physics (e.g. the self-avoiding walks), fractal 
patterns playing important roles in non-equilibrium statistical mechanics models (e.g. the 
loop-erased random walks and the uniform spanning trees for sandpile models and forest 



fire models showing self-organized criticality), and so on (4-6]. The theory is based on two 
branches of mathematics, the complex function theory [3] and the stochastic analysis {2], 3| 
and is strongly connected with the conformal field theory js-10]. 

As well as by wideness of applications and by richness of mathematics, we are attracted 
by simple setting of the theory ; Consider a complex plane C. (i) First we consider a motion 
of Brownian particle on the real axis M. We assume that it starts from the origin at time 
t = and the diffusion constant is given by k > 0. If we write the position of the diffusion 
particle on R at time t > as Ut, then (Ut) = and (U?) = Kb for t > 0. Usually we 
denote the position of a one-dimensional standard Brownian motion (BM) at time t by B t , 
for which (B t ) = B Q = 0, (B 2 ) = t,t > 0. The BM has the scaling property such that for 
any constant c > 0, the distribution of the position of BM at time c 2 t is equal to that of the 
position of BM at time t multiplied by a factor c, that is, the equality B c 2 t = cB t holds in 
distribution. Then we can give the above U t by 

U t = VnB tl t > 0. (1) 

(ii) Then we solve the following partial differential equation for a complex function ft(z) on 
the upper half complex-plane HI = {z G C : Im z > 0}, 

dMz)_ = _df 1 (z^J2_ 
dt dz z-U t U 

under the initial condition fo(z) = z. (iii) Note that the boundary of EI consists of the real 
axis K and an infinity point. When z G HI approaches the special point Ut on R, the position 
of the BM, the RHS of ()2]) diverges. If we trace the image of this singular point 

lt = ft(U t ), t>0 (3) 



we will have a curve 7(0, t] = {7^ : < s < t} in EI starting from the origin. For each s > 
the curve 7(0, s] and the region enclosed by parts of the curve and the real axis R should be 
eliminated from H in order to continue to solve Eq. (f2|) for t > s. 

For any deterministic simple curve 7(0, t] in H, t > 0, Loewner proved that there exists a 
real-valued function U t and a conformal map f t , which is one-to-one from EI to EI \ 7(0, t] = 
{z G H : z 7^ 7 S , < s < t}, the upper half complex-plane with a slit 7(0, t], such that f t 
and U t solve Eq.© with the condition (13|). The equation (]2]) is called the Loewner equation 
[3, Schramm considered an inverse problem: given Ut on R and derive a curve 7(0, t] by 
solving Eq.([2]). Since he gave U t by a BM as Eq.flTJ), ^ t given by Eq.(E]) performs a stochastic 
motion and the obtained curve 7(0, t] is statistically distributed in HI [lj. Equation (T2J) 
driven by a BM with variance nt is called the stochastic Loewner equation or the Schramm- 
Loewner evolution (SLE) 12|, and a random curve 7(0, t], t > is called the SLE curve with 
parameter k (the SLE K curve). 

Although the change of diffusion constant k causes only quantitative change of the driving 
function U t , that is, scale/time change \JHB t = B Kt in distribution, it does qualitative change 
of SLE K curves. When < k < 4, the curve is simple (with no self-intersection) with 
7(0, 00) C EI. When 4 < k < 8, the curve is self-intersecting and 7(0, 00) hits the real axis 
R infinitely many times, but it is not dense on EI : 7(0, 00) D H 7^ H with lim^oo |7 t | = 00. 
And when k > 8, it will cover whole of H in t — > 00 J2j3|. The fractal dimension (Hausdorff 
dimension) of the SLE K curve is determined as 13|, llij 



d K 



1 + k/8, 0<k< 
2, k>8. 



(4) 



Moreover, if k is chosen to be a special value, SLE K curves provide the statistical ensembles of 
continuous limits of random discrete paths studied in statistical mechanics models exhibiting 



critical phenomena or in fracta 
critical percolation model 15|, 



models on lattices. For example, the value k = 6 is for the 



161 ] . Effective methods for numerical simulations of SLE K 



curves by computers have been reported [17H19|. 

Today we can learn from mathematics literatures that the Loewner equation (T5]) for a 
deterministic function Ut had played important roles in the complex function theory even 
before Schramm introduced its stochastic version [7]. We shall say, however, that this 
equation has not been familiar to us, physicists. More familiar differential equation to us in 
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the complex calculus is the one, whose solution gives a conformal transformation from EI to 
the interior of a polygon on the complex plane with mapping the real axis M to a piecewise 
linear boundary of the polygon, called the Schwarz-Christoffel transformation (SCT) (see, 



for example, 
case of SCT. 



20(|). So here we try to discuss the Loewner equation (J2J) by using a special 




-x L x° x R 



FIG. 1: SCT generating a straight slit in H 

Let < a < 1. Consider a conformal map / from H to the upper half complex-plane 
with a straight slit starting from the origin, H \ {a slit}, where the angle between the slit 
and the positive direction of the real axis is an as shown in FigJTJ Since EI with the straight 
slit can be regarded as a polygon with the interior angles (1 — a)n on the left side of the 
origin, 2n around the tip of the slit, and an on the right side of the origin, for any length of 
a slit, the conformal map is given as an SCT, which solves the differential equation 



df(z 



= b 1 (z + x L y a (z-x°)(z-x R ) a -\ (5) 



dz 

where b\ is a complex number and real numbers satisfying the inequalities, 

x L > 0, x R > 0, — x L < x° < x R . By this transformation, both of the points — x L and x R on 
R are mapped to the origin, /(— x h ) = f{x R ) = 0, and x° to the tip of the slit. See FigJTJ 
We have found the general solution of (J5j) expressed by 

f(z) = b + b 1 (z + x h f- a (z - x R ) a 

C x° -x R + a(x L + x R )T(l-a) ^ f x L + z \) 

>Ml + i s ^^7 {F M, 2 -a;- T ? H, 6 

\ x L + x R r(2 - a) V x + x J J 

where b is a complex number, T(z) = J °° e~ u u z ~ x du (the gamma function), and F(a, (3, 7; z) 
is Gauss's hypergeometric function F(a, z) = T (j) / {T (a)T (f3)} Y^=o^( a + n )^(.P + 
n)z n /{r(7 + n)n\}. We impose the hydrodynamic condition 



z 



— > 1 in z — > 00. (7) 



Then bo = 0,b\ = 1 and x° — x R + a(x L + x R ) = 0, and we have f(z) = (z + 
Eq.® is rewritten in this 

df(z) 2 2f(z) 



x \ J \l-ot( z _ x K\<x. 



(8) 



dz z — x° (z + x L )(z — x R ) 
We then introduce a parameter "time" t > and assume that x L and x R , and thus also a; , 
depend on t by setting 

f t (z) = (z + x);) i -«(z-xfr. (9) 

The differential of ft with respect to t is written as 

df t (z) 2A(z,t)f t (z) 
dt (z + x\)(z — xf) 

withA(z,t) = [-(l-a){z-x R )dx\ / dt + a{z + x\)dx R / dt}/2. Let = 2ct> 3 ,x R = 2^/cwith 
constants c > and (3. Then we can see that, if and only if c = yJa/{\ — a) and f3 = 1/2, 
A(z,t) becomes independent both of z and t; A(z,t) = 1. Combining this observation with 
Eq.® gives the following result: For < a < 1, the SCT 



l-a 



^> = HvrVv (W^J (10) 

is not only a solution of the Schwarz-Christoffel differential equation (J5J), but also of the 
Loewner equation (J2]) with the driving function 

U? = x° t = xf - a{x)> + x R ) 
[kFI if a < 1/2 



_ (11) 
k q £ if a > 1/2, 



where 



4(1-2^ 

a l-a v 7 



As time t goes, the straight slit performs as an "evolutionary boundary" of the image of HI 
by f t , in which the tip of the slit fl3]) is evolving as 

it = mun 

= 2 ( — -^J e^Vi, t > 0. (13) 
One can observe the equality 



u? = ±J((V^B t y), t>o, (14) 



since (Bf) = t. That is, the driving function Uj* for the above SCT (TTTT) is the positive or 
the negative root square of the squared average of the random driving function (JT|) of the 
SLE. So if we are able to introduce fluctuations into the SCT systematically, we can draw 
approximations of SLE curves on H. It may be the basic idea to simulate SLE curves by 
dividing a time period (0,t] into n small intervals {(tj-i,tj] : j = l,2,...,n} by setting 
= to < t\ < ■ ■ ■ < t n = t and the above single SCT is replaced by an n-multiplicative map 



of infinitesimal SCTs with sufficiently large n 
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On the other hand, we know the fact that the diffusion property of BM can be observed 
in long-time asymptotic behavior of a simple discrete-time stochastic process, random walk 
(RW). Let <7j, j = 1, 2, 3, . . . be independent and identically distributed (i.i.d.) random vari- 
ables taking values <jj = 1 with probability 1/2 and Oj — — 1 with probability 1/2. Consider 
a simple symmetric RW on the one-dimensional lattice Z = {. . . , —2, — 1, 0, 1, 2, . . . } start- 
ing from the origin at time n = 0. We denote the position of the random walker at time 
n = 0, 1, 2, . . . by w n . Then, wo = and 

n 

w n = ^2<Tj, n= 1,2,3,.... (15) 

In the present paper, we consider an SCT as a functional of a random variable a and consider 
an Iterative system of SCTs (ISCTs) driven by RW. By this system, each time series of steps 
(<7i, (72, ... ) of RW is mapped to a series of points £ 2 > • • • ) in HL Let W n and S n , n > 
be the interpolations by straight lines of w n and £ n , n = 0, 1, 2, . . . , respectively. Since for 
T > 0, {W^i/v 7 ^ : < t < T} converges to {B t : < t < T} as iV oo in distribution, 



C (0, T] = ■{ —j=i < t < T , T>0, (16) 



which we call the rescaled ISCT path with up to time T, will converge to an SLE curve 
up to time T, 7(0, T], in distribution as iV — > 00. Figure [2] shows the rescaled ISCT paths 
^(0, 1] = {Z m /y/N,0 < t < 1} with iV = 5 x 10 4 for k = 2 and 6, which are drawn 
by interpolating the series of points {£j/\N}jLo by lines. They seems to approximate the 
SLE K curves 7(0, 1] very well. In the present paper, we will show that, for sufficiently large 
N, nas the same statistical properties as the SLE curve 74 has. 

The paper is organized as follows. In Sec. II, we define the ISCT driven by random walk 
with a parameter k = n a defined by (jT2l) and the ISCT K paths. In Sec. Ill, by observing 
the behavior of ISCT paths S(0,n] for small values of n, we show that, even if we use the 




FIG. 2: Approximations of the SLE K curves up to time T = 1 expressed by the ISCT K paths for 
k = 2 and k = 6. 

same realization of RW as a driving function, the ISCT K paths with large values of k exhibit 
much more complicated motion on H than those with smaller values of k. In Sec. IV, we 
report the properties of the rescaled ISCT paths ( N (0,T] with T = 1 based on large scaled 
computer simulations with iV ~ 10 4 ~ 10 5 . We evaluate the fractal dimensions dj SCT of 
the iV — > oo limits of C N {®i 1] f° r several values of k, and the dependence of c?f SCT on k 
is compared with that of the fractal dimensions d K of the SLE curves given by (jl]). We 
also show that, for 4 < k < 8, the generalized version of Cardy's formula of SLE K curves 
[2I. I21I] will be applicable to the rescaled ISCT K paths ( N (0,T\, if and T are sufficiently 
large. Section V is devoted to give concluding remarks. Appendix A is prepared for giving 
recurrence relations, which will be useful to analyze the ISCT paths. 



II. ISCT DRIVEN BY RW 



Noting that Eq. (}T2]) is solved for a as a = [1 ± ^/k/(k + 16)]/2, we set 



a en 



1-a 



K+1Q_ 

for k > 0, cr G {—1, 1}, and define the SCT as a functional of a random variable a by 



1— a K (c) / i ^ a K (a) 



(17) 




^ * (z-2J 1 - a ;W 1 . (18) 



1 — a K {a) J \ V oi K [a) 

Given one step <j\ of RW on Z, we consider an SCT, F£ (z), which is a conformal map 
from HI to the region HI with a straight slit. The straight slit starts from the origin and ends 
at the tip located at 



6 = i^(v^i) 



1 -a K (ai) 



l/2-« K (<Tl) 



3 ia K (cri)7r 



(19) 

Next assume that two steps of RW, (ci, cr 2 ), is given. We transform HI by an SCT, F£ 2 . The 
image of HI, F* (H), is the region HI with the straight slit, which starts from the origin and 
ends at £2 = -^(V^ " 2 )- Then we consider the transformation of the region HI with this 
straight slit by another SCT, i 7 ^. By this SCT, a straight slit from the origin to the point 
£1 is generated as shown by (fl9|) . The image by F£ of the straight slit between the origin 
and £2 m i s 5 however, no longer a straight line but a curvy one. It starts from £1 and 
ends at 

£ 2 = F; i ( v ^a 1 + £«) 

= F«(v^i + F*(v^T a )). (20) 

We have then a set of two points (£1, £2) in HI. Set n > 1 and now we assume that a realization 
of RW on Z up to time n is specified by the series of n steps, <r(n) = (01, 02, . . . , cr n ). Let 

= J?(vfo7 + *) (21) 
for k > 0, cr G { — 1, 1} and z G HI. We perform the following iteration of SCTs, 

= 5- (5«(.-.(^ n (z)) •••))• (22) 



Then we have a curve consisting of a straight slit between the origin £o = and £1 and n — 1 
segments of curvy slits, which are sequentially connected at £j, 1 < j < n — 1, and the tip 
is at £ n . See Fig. [3J In the present paper we call Sg., n Jz) the iterative Schwarz-Christoffel 
transformation (ISCT for short) driven by RW specified by <r(n) = (a±, . . . , a n ). 

*&) 

f S-X 




FIG. 3: Iteration of SCTs 

For n — 1, 2, 3, . . . , we define the points in H by 

£» = S&(„)(0). (23) 

We set £o = 0. The sequence of points (£o,£i,£2, • • • , £n) is interpolated by straight lines. 
We call it an ISCT path in EI and denote it by H(0, n}. In other words, each realization <x(n) 
of RW on Z is mapped to a path 2(0, n] on H by the ISCT. 

For a given n > 1, we introduce the following recurrence relations for a series 

(AO) A 1 ) An)\ 
yc,n ,Zn j • • • j Zn )i 

e +1) = S: n _M k) )> fc = 0,l,2,...,n-l (24) 
with = 0. Since = 1,2,3, ...,n are i.i.d., the recurrence formula (1241 is useful 
to calculate the position £ n , which is given by Moreover, if we introduce a parameter 
8 G (0, 7r/2) with the relation 



cos6=J K (25) 



K+16' 



171) is written as 



a K (cr) = -(1 - crcosi 



sin 2 (#/2), a = l 

(26) 

cos 2 (fl/2), a = -1, 



and (fT9|) is given by 

6 = 5^(0) 

2(cot(fl/2)) cos2 V 7rsin2 ( / 2 \ <n = 1 
2(cot(#/2)) cos2 V 7rcos2 W 2 ), a 1 = -l. 

Note that £i(er = —1) = — (£i(er = 1))*, where * indicates complex conjugate. The recurrence 
relation f)24p is then written as 

, n\ cos 2 (0/2) , a\ sin 2 (6/2) 

Cr^^ + ^cot-j ^i fc) -a n _ fc tan-J , k = 0, 1, 2, . . . , n - 1 

(28) 

for cr ra _fc = ±1 with ^i '' = 0. The expressions for the relations between two real components 
of complex variable ^n +1 ^ and those of £n are given in Appendix A. 



III. NETWORKS ON M GENERATED BY RW PATHS 

For a fixed time period n > 0, consider a collection of all realization of steps of RW on Z, 

JT( n ) = | CT ( n ) = a 2 , . . . , 0-„) : a, E {-1, 1}, 1 < j < n}. (29) 

The total number of realizations of RWs is 1^(^)1 = 2 n . We note that each realization of 
RW up to time n > is represented by a directed path on a squared lattice in a triangular 
region in the spatio-temporal plane, 

A n = {(j,k) G Z X {0, 1,2, ...} : < k < n,-k < j < k, j + k = evenj. (30) 

Here any edge connecting nearest neighbor vertices in A n is assumed to be directed in the 
positive direction of the time axis; (j, k) — > (j ± 1, k + 1), and each path is a sequence of 
edges all in the positive direction, starting from (0,0) to (j, n) with — n < j < n. Figure H] 
shows A4 and an example of a realization of RW with <r(4) = (1, 1, —1, 1). 
By collecting all ISCT paths up to time n, S(0,n], we have a network on H, 

AC = {S(0,n] : = S£ (i) (0),l < j < n,a(n) G JT(n)}. (31) 

Figure [5] shows N% for k = 2, 4, 6 and 8 for n = 4. As the parameter k increases, the network 
becomes spreading wider in H. There the ISCT path corresponding to the realization of RW 
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space 



FIG. 4: Lattice A4 and one realization of RW with tr(4) = (1, 1, —1, 1). 

cr(4) = (1, 1, —1, 1) shown in Fig. H]is indicated by a bold line for each value of n. In Fig. 
[6J we compare the ISCT paths with (b) k = 2 and (c) k — 8 both obtained from the same 
realization (a) of RW with n = 20 steps, where the networks up to n = 10, Afi , are also 
shown in the background for each k. We can see that the ISCT path with k = 8 is much 
more complicated than the path with k = 2. As shown by Figj5] and Figj6j the network 
is bounded by the rightmost path H max (0, n] generated by <x max (n) = (1, 1, . . . , 1) and 
the leftmost path S min (0, n] generated by cr min (n) = (—1, —1, . . . , —1). The height of ~™ ax 
and S™ in , H£ = ImS™ ax = ImS™ in , is observed to converge to a positive constant in 
n — > 00. The numerical values are given by = 0.014(k = 1), 0.010 (k = 2), 0.0084 
(k = 3), 0.0074 (k = 4), 0.0067 {k = 5), 0.0062 (k = 6), 0.0058 {k = 7), and 0.0055 {k = 8). 

IV. NUMERICAL ANALYSIS OF CURVES GENERATED BY ISCT 

For a given number of steps n = NT of RW, we have defined the rescaled ISCT path 
( N (0,T] by Eq. ljTB"]) . In particular, ^(0, 1] is obtained by interpolating the series of points 
{£j/\ /r N}f=o by straight lines. We have studied statistical properties of the rescaled ISCT 
paths based on the numerical data of large scaled computer simulations. 



11 




FIG. 5: Networks Af* on EI for n = 4 with k = 2, 4, 6, 8. The ISCT paths corresponding to the 
same realization of RW <r(4) = (1,1,— 1,1) are indicated by bold lines. 
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FIG. 6: (a) A realization of RW with time step n = 20. (b) The ISCT path for n = 2 generated 
from the RW shown in (a) with the network Nf - (c) The ISCT path for k = 8 generated from the 
same RW shown in (a) with the network Mf . 
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A. Fractal Dimensions 



Figure [7] shows a log-log plot of the box counting of segments of ( N (0, 1] with respect to 
the box sizes for k = 4 with N = 5 x 10 4 . As shown by this figure, the data for any k can 
be fitted by a straight line very well and we can evaluate the approximate values of fractal 
dimensions for finite N. 

10 
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FIG. 7: Log-log plot of the box-counting data for the restricted ISCT paths for n = 4. 

The evaluated values up to at most N = 3 x 10 5 are then plotted versus 1/N in FigJSJ 
For each evaluation we prepared M = 20 samples, and the ranges of scattering of results are 
shown by error bars in the figure. There N — > oo limits are extrapolated by three-parameter 
fittings ; d = a + a±/N + a 2 /N 2 . The obtained values ao by the extrapolation are denoted 

by ^fscT- 

Figure M shows dependence of df SCT on k. The Hausdorff dimensions of the SLE K curves 
given by (jl]) are also shown by a dotted line. Systematic deviation is found between d K for 
the SLE curves and df SCT numerically evaluated for the rescaled ISCT paths. We observe 
in FigJH] that the approximate value of fractal dimension for finite N is increasing as N 
is increasing, and the ratio of increment becomes larger as k approaches the value 8. So 
we expect that the deviation will be systematically reduced, if we can perform numerical 
simulation for larger TV's and make appropriate extrapolation to the N — > oo limit. 
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FIG. 8: Extrapolation of fractal dimensions in the N — > oo limit. 




FIG. 9: Dependence of the fractal dimensions df SCT of rescaled ISCT paths is shown. The Hausdorff 
dimensions of the SLE K curves given by Eq.(j3]) are also plotted by a dotted line. 

B. Generalized Cardy's Formula 

Here we first consider the SLE curve 7(0, 00) = '■ < t < 00} in the case 

4 < k < 8. (32) 



In this case, jt starting from the origin will hit the real axis K infinitely many times [2 
For x > we can define 



3]. 



t*(x) = the first time, when 74 hits a point in [x, 00) on R. (33) 
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Then, 7^( x ) is the leftmost point in the interval [x, oo), at which jt hits R. Note that, if 
k < 4, t*(ar) = oo, and if k > 8, 7t»( x ) = x, with probability one. For (132]) . 7^(2;) has a 
nontrivial distribution. For each e > 0, we observe whether ^(x) < x + £ or £t„,( x ) > x + e. 
Figure [TU1 illustrates the former case. For the SLE K curves with (|32[) . the following formula 
is established. (See Proposition 6.34 in [2|.) 

t(4/k) r £/(x+£) 

W) Jo 



P(7t.(x) < x + e) 



r(8/«- i)r(i 



u °/«-2(l 



r(4/«) 



8/k-1 



e 



I -,- - 1,- 

Kj Kj Kj X ~\~ s 



where T(z) is the gamma function and F(a, ft, 7; z) is Gauss's hypergeometric function. This 
formula can be regarded as a generalization of Cardy's formula, since the original formula 



corresponding to f l34j) with k = 6 was derived by Cardy 



21] for the "percolation exploration 



process" in the critical percolation model, and then the continuum limit of that process was 
proved to be described by the SLE curve with k = 6 [lsl, Q. For < e <C 1, the above 



formula gives a power law 
with the exponent 



S(k) 



— K 



K 



(35) 
(36) 








x 



| X + s 
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FIG. 10: Illustration of the point r Yt»(x) on the rea ^ ax i s ^ f° r an SLE K curve with 4 < k < 8. 

Now we consider the ISCT paths. For integers N and T with N ^> 1, prepare a realization 
of RW represented by cr(NT) = (cri, cr 2 , . . . , cjnt)- For each n = 1,2,..., NT, by using the 
data (<Ti, cr 2 , . . . , <j n ), we calculate the position in EI following the recurrence formula (1241) . 
As noted at the end of Sec. Ill, Im£ n > for any n > 1. So we set a small value h > and 
look for the event 

:N / Im£ n Re£ n 



EFJn) = < < h and — ^ > x } . (37) 
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We define 

rih(x) = min |n : 1 < n < NT, S^ x (n) occurs j. (38) 

If n h (x) < NT, we define t h (x) = n h (x)/N and calculate the value Re^, x -, = 
ReS^^/V^/V = Ke^ h ^/y/N. If n h (x) > NT, that is, the event ( 157)) does not occur 
for the given cr(NT), then th{x) > T. The probability distribution function for the rescaled 
ISCT paths, which corresponds to P(jt*( x ) < x + e), may be given by 

lim lim lim p(lle(? (x) <x + e, t h (x) < t) = P(Cux) <x + e), (39) 
where C(0, oo) = Xxm^^Xvax^^ C^ {^,T\. 

1 i ■ i 




OO & 1 O001 0.001 0.01 0,1 1 



epsilon/(x+epsilon) 



FIG. 11: Log-log plots of the numerical evaluations of the probability P(ReQ^ a ,x < x+e, th(x) < 1) 
versus e/{x + e) for the rescaled ISCT paths. 

In numerical calculations, we have set iV = 10 4 and T = 1 and prepared M = 1000 
realizations of RW, which are randomly generated. Since T = 1, we are allowed to consider 
only small values of x and e. The threshold value h is wanted to be small, but for finiteness 
of N, it should be positive. In FigJTTj we show log-log plots of the numerical evaluations of 
P{Re( t N h(x) <x + e, t h (x) < 1) for « = 5 with x = 0.8, h = 0.1, « = 6 with x = 0.5, h = 0.1, 
and k = 7 with x = 0.8, h = 0.2. For finiteness of N and smallness of the number of samples 
M, data scatter for small e. We find, however, power-law behaviors in the intermediate 
regions of e; 

P(Cu(x)<x + e)^e s . (40) 

By linear fitting as shown in FigJTT) we have evaluated the values of exponent 5. The results 
are plotted in FigJT21 where Eq. (l36)) derived from the generalized Cardy's formula is also 
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shown by a curve. The good agreement implies that in the proper limit (139|) will also follow 
the generalized Cardy's formula in the parameter region (|32|) . 




kappa 



FIG. 12: Numerical evaluations of the exponent 5 in the power law (|40p for k = 5,6, and 7. The 
curve shows Eq. (j36p of the generalized Cardy's formula. 



V. CONCLUDING REMARKS 

In the present paper, we have proposed an algorithm, which generates a random discrete 
path 2(0, NT] on the upper half complex-plane EI as a functional of a path of random walk 
W(0, NT] on the real axis R for integers N and T. The system has one parameter k > 
and we call the path an ISCT re path. We have studied the rescaled ISCT K path defined by 
( N (0,T] = H(0, NT\/y/~N for large N by computer simulations. The numerical analysis of 
the distributions of ( N (0,T] supports our expectation that the limits of the rescaled ISCT K 
paths C(0, oo) = limT-).oo hmjv^oo ^(0, T] will have the same statistical properties as the 
SLE K curves have. 

The rescaled ISCT K paths ( N (0,T] can be regarded as discrete approximations of the 
SLE K curves. In this sense, the present study could be included by the previous numerical 



work 17j,ll8]. In this paper, however, we have emphasized on our interest in the ISCT itself as 
a simple algorithm to generate complicated discrete dynamics of a point on H. Dependence 
on the parameter k of complexity of the SLE curves is demonstrated by dependence of 
complexity of the network of the ISCT paths on the angle om of the slit generated by a 
single SCT. 
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We have learned that stochastic analysis is necessary and useful to study statistics and 
stochastics of the SLE K curves {2, 3|. Although we have reported only numerical study in 
this paper, we hope that the combinatorics and statistical mechanics methods developed for 
solvable models on lattices will be useful to analyze statistics and stochastics of the ISCT K 
paths on H, since they are functionals of simple random walks in one dimension. 
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Appendix A: Recurrence Relations 

Let £ n k) = 2re l(f> , r > 0, < < vr for k < n. Then (|28]) gives 

^ fc+1) = 2R^ k exp(i* ffB _J, a n . k = ±1 (Al) 



with 



R a = 2 ( r 2 + 2ar cos <fi cot - + cot 2 - 



g (cos 2 (fl/2))/2 



q q\ (sin 2 (0/2))/2 



x yr 2 - 2ar cos tan- + tan 2 . (A2) 

/ rcos0 + <7-cot(0/2) \ 2 9 

$ CT = arccos — . cos - 

\y/r 2 + 2arcos0cot(fl/2) + cot 2 (fl/2) J 2 

(r cos d> — a tan (9/2) \ . 9 9 , . . 

v ' ' = sm 2 - (A3) 
y/ r 2 -2<xrcos0tan(fl/2) + tan 2 (fl/2) J 2 

When we set ffi = 2{x + iy), x G M, y > 0, k < n, the above gives 

= 2(X CTn „ fc + iY an _ k ), a n . k = ±1 (A4) 
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with 
X„ -- 



q q\ (cos 2 (0/2))/2 

x 2 + y 2 + 2ax cot — h cot 2 
y 2 2 



2 2 ^ 2 ^ 

x + y — 2ox tan — h tan - 
y 2 2 



(sin 2 (0/2))/2 



X COS 



x + acot (9/2) 



2$ 

arccos | — ' ; ' = | cos - 

^x 2 + y 2 + 2axcot (9/2) + cot 2 (9/2) I 2 



+ arccos 



x — crtan (9/2) 



'x 2 + y 2 - 2(xxtan (9/2) + tan 2 (9/2) 



x 2 + y 2 + 2(jx cot - + cot 2 - 







( cos 2 (e/2))/2 



• 2° 

sin - 
2 







(A5) 



x 2 + w 2 — 2crx tan - + tan 2 - 
2 2 



(sin 2 (9/2))/2 



x sin 



arccos 



x + crcot (6>/2) 







cos 



V^ 2 + y 2 + 2as cot (9/2) + cot 2 (0/2) / 2 



+ arccos 



.t - crtan (9/2) 



sin 



v/x 2 + y 2 - 2crx tan (0/2) + tan 2 (9/2) J 2 



(A6) 
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